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HIGHLIGHTS 


•  Electrokinetic  energy  conversion  is  investigated  experimentally. 

•  High  electrokinetic  energy  conversion  efficiency  in  Nation  membrane  is  found. 

•  Power  density  and  efficiency  increases  with  temperature. 

•  Temperature  dependence  of  pore  size  in  Nation  are  estimated  from  transport  properties. 
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In  the  present  study  we  investigate  the  transport  properties  of  Nafion®  117  membranes  in  temperatures 
ranging  from  ambient  temperature  up  to  70  °C.  The  hydraulic  permeability,  streaming  potential  and  ion 
conductivity  have  been  measured  as  function  of  temperature  in  0.03  M  LiCl  solutions  in  purposely 
designed,  non-conductive  set-ups.  In  particular,  the  apparent  activation  energies  of  the  processes  have 
been  retrieved:  29.4  kj  mol-1,  9.3  kj  mol-1  and  22.9  kj  mol-1  for  the  hydraulic  permeability,  streaming 
potential  coefficient  and  ion  conductivity  respectively.  Based  on  the  knowledge  of  the  temperature 
dependence  of  these  three  independent  properties  the  figure-of-merit  of  the  electrokinetic  energy  con¬ 
version  process  has  been  calculated  obtaining  a  monotonous  increase  of  the  efficiency  with  temperature. 
At  70  °C  the  electrokinetic  efficiency  is  rather  high  about  26.6%:50%  higher  with  respect  to  the  one  found 
at  room  temperature.  The  electrokinetic  transport  properties  were  also  used  to  esteem  the  average  pore 
size  of  the  water  channels  in  the  polymer  matrix  resulting  in  pore  diameters  ranging  approximately  from 
2.0  (25  °C)  to  2.8  nm  (70  °C). 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  increasing  demand  of  energy  and  the  need  of  decentralized 
power  generation  systems  combined  together  with  the  urgency  of 
emission  reductions  encouraged  the  researchers’  interests  in 
developing  new,  efficient  and  more  environmental  friendly  energy 
production  systems.  Polymer  electrolyte  fuel  cells,  PEMFCs,  are 
potentially  among  the  most  efficient  devices  for  power  generation, 
and  a  large  amount  of  specific  literature  has  been  published  in  the 
last  decades  focusing,  especially,  in  developing  new  solid  polymer 
electrolyte  materials  with  enhanced  properties  [1—6].  However, 
despite  the  efforts  devoted  in  this  direction  the  reference  mem¬ 
brane  for  fuel  cell  applications  still  remains  Nafion®,  produced  by 
DuPont,  due  essentially  to  the  high  chemical,  mechanical  and 
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thermal  stability  and  high  proton  conductivity  during  fuel  cells 
operation.  Another  attractive,  still  not  well  investigated,  way  of 
energy  conversion  in  ion  conductive  membranes  is  the  direct 
conversion  of  kinetic  energy  into  electrochemical  energy.  The 
electrokinetic  energy  conversion  has  been  studied  both  theoreti¬ 
cally  7-10]  and  experimentally  focusing  the  attention  mostly  on 
the  properties  of  straight  nanochannels  with  reported  efficiencies 
up  to  10%  [11-16].  In  literature,  a  few  studies  of  the  electrokinetic 
properties  of  ion  conductive  membranes  have  been  reported.  A 
non-mechanical  pressurization  system  based  on  Nafion®  117,  in  the 
following  referred  to  as  Nafion,  in  tetrapropylammonium  iodide 
and  iodine  in  DMF  solutions  has  been  proposed  by  Evans  et  al.  [17], 
while  we  have  reported  the  electrokinetic  properties  of  Nafion  at 
room  temperature  in  a  previous  work  [18].  However  a  more  sys¬ 
tematic  experimental  investigation  of  the  electrokinetic  properties 
of  ion  conductive  membranes  is  still  missing.  The  basic  mechanism 
behind  the  electrokinetic  power  generation  is  the  potential 
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difference  generated  by  a  pressure  gradient  applied  across  a 
membrane,  i.e.  streaming  potential,  separating  two  compartments 
filled  with  aqueous  solutions  having  the  same  salt  concentration. 
The  streaming  potential  is  a  consequence  of  the  coupling  between 
ions  and  solvent  (hydration)  and  whenever  a  pressure  difference 
drives  a  volume  flux  through  the  membrane,  ions  are  dragged 
along,  whereby  creating  a  potential  difference.  One  of  the  key  ad¬ 
vantages  of  electrokinetic  energy  conversion  is  that  moving  me¬ 
chanical  stages  are  not  required.  Nonetheless,  high  conversion 
efficiency  along  with  high  power  density  are  fundamental  param¬ 
eters  in  comparisons  to  e.g.  electromagnetic  actuators  that  repre¬ 
sent  the  state-of-the-art  of  commercially  available  units  that 
convert  kinetic  energy  into  electrical  energy. 

Based  on  measurements  of  the  hydraulic  permeability  (kh), 
streaming  potential  coefficient  (v)  and  ion  conductivity  {a)  the 
electrokinetic  figure-of-merit  /?ek  =  v2  •  (t/k h  is  retrieved  and  can  be 
used  for  calculation  of  the  electrokinetic  conversion  efficiency 
1/EK  =  [(l+ftx)1/2  -!]/[(!+  ftac)1/2  +  l].  For  a  more  detailed 
discussion  the  reader  is  referred  to  refs.  [7,8,16,19]. 

In  previous  studies  we  demonstrated  that  the  maximum  effi¬ 
ciency  of  18  ±  2%  for  Nation  [18]  is  significantly  higher  with  respect 
to  the  non-  or  low-charged  membranes  [16,18],  for  which  a 
maximum  of  about  7%  for  Cyclopore  and  Nucleopore  membranes 
has  been  found  16].  These  results  indicate  that  ion  conductive 
membranes  could  represent  a  possible  pathway  for  high-efficiency 
electrokinetic  energy  conversion.  The  high  electrokinetic  efficiency 
observed  in  ion  conductive  membranes  with  respect  to  non- 
charged  membranes  could  be  explained  by  the  high  surface 
charge  density,  aq  ~  0.2—1  C  m-2  for  Nation  [20,21],  of  the  nano- 
porous  hydrophilic  network.  In  its  swollen  state,  an  ion  exchange 
membrane  can  be  depicted  as  a  hydrophobic  PTFE-like  matrix  and  a 
network  of  more  or  less  interconnected  water  nano-channels  in 
which  the  transport  of  water  and  ions  takes  place  [22-27].  Besides 
the  fixed  charge  density,  the  electrokinetic  properties  of  such  a 
system  strongly  depend  on  the  average  pore  dimensions  and 
channel  tortuosity  [28,29].  Several  ways  to  modify  the  pore  struc¬ 
ture  of  Nation,  to  increase  the  ion  exchange  capacity,  the  water 
retention  and  the  mechanical  and  thermal  stability  have  been 
proposed  in  the  literature.  This  includes  e.g.  addition  of  hygroscopic 
materials  (AI2O3,  Si02,  clays,  silicophosphate  gels  etc.)  obtaining 
inorganic-polymer  composite  membranes  [30],  use  of  dendrimers 
[31  ] .  On  the  other  hand  it  is  well  known  from  the  literature  that  by 
increasing  the  temperature,  Nation  membranes  show  higher 
dimensional  swelling  [32],  liquid  water  uptake  [33-35]  and  hy¬ 
draulic  permeability  36,37].  These  latter  experimental  evidences 
indicate  a  more  open  structure  with  higher  water  content  and 
larger  average  pore  size.  It  is  anticipated  that  all  the  three  transport 
properties  that  enter  the  expression  for  (3ek  =  v2  •  oV^h  increase  with 
the  pore  dimensions  and  thus  with  temperature.  Still,  the  expected 
temperature  dependence  of  /?EK  is  not  known  and  depends  on 
whether  the  nominator  (y2-o)  has  a  larger  increase  with  tempera¬ 
ture  with  respect  to  the  denominator  (kh).  Therefore  the  main 
objective  of  this  work  is  to  experimentally  determine  the  temper¬ 
ature  dependence  of  kh,  v  and  a  in  Nation  membranes  and  to  get 
more  detailed  insight  on  the  electrokinetic  conversion  efficiency  as 
function  of  temperature. 

2.  Materials  and  methods 

2.1.  Membranes 

Nation  117,  with  a  nominal  thickness  of  178  pm,  was  bought 
from  Ion  Power  Inc,  US.  Before  use  membrane  samples  were  pre¬ 
treated  by  boiling  in  demineralized  water.  The  membranes  were 
then  heated  in  3wt%  of  H2O2  water  solution  for  60  min  at  T  ~  90  °C, 


transferred  in  boiling  demineralized  water,  and  then  heated  up  to 
90  °C  for  60  min  in  0.5  M  H2SO4  water  solution.  Finally  the  mem¬ 
branes  were  washed  in  warm  ( ~  70  °C)  demineralized  water  for  ten 
minutes;  the  last  step  was  repeated  four  times.  After  pretreatment 
the  membrane  samples  were  stored  in  0.5  M  LiCl  solutions  to 
obtain  Nation  in  Li+  form.  All  measurements  have  been  carried  out 
in  0.03  M  LiCl  aqueous  solutions,  and  to  avoid  any  increase  in  the 
concentration  in  the  measurement  cells  the  membranes  were 
carefully  rinsed  in  0.03  M  LiCl  before  being  tested.  The  membrane 
thickness  in  both  the  dry  and  swollen  states  was  measured  with  a 
digital  micrometer  (Diesella,  resolution  1  pm). 

2.2.  Hydraulic  permeability  and  streaming  potential  apparatus 

The  experimental  set-up  for  the  measurement  of  the  hydraulic 
permeability  and  streaming  potential  as  a  function  of  the  temper¬ 
ature  is  schematically  shown  in  Fig.  1.  It  consists  of  two  jacked 
borosilicate  glass  chambers  of  about  400  cm3  volume  separated  by 
two  membrane  holders  made  in  Teflon.  The  custom-made  holders 
had  4  ports  each  for  instruments’  connection  (e.g.  thermocouples 
and  electrodes).  The  temperatures  in  the  two  cells  were  indepen¬ 
dently  controlled  by  two  different  water  baths  (Eco  Gold  Lauda  DE, 
temperature  stability  0.01  °C)  feeding  the  chambers’  jackets,  and 
monitored  through  two  thermocouples  type  I<  (Omega  US,  1/ 
16  inch)  placed  at  approximately  5  mm  from  the  membrane  sur¬ 
faces.  To  reduce  concentration  polarization  at  the  membrane/so¬ 
lution  interface  and  to  have  a  more  uniform  and  constant 
temperature,  two  inert  magnetically  coupled  stirrers  (Biichi  CFI, 
Cyclone  i)  have  been  employed.  The  stirrers,  having  all  the  wet 
parts  in  PFA  and  glass,  are  rated  to  work  up  to  6  bar  of  pressure  and 
temperature  up  to  180  °C  and  they  were  connected  through  GL  32 
ports  to  the  two  Teflon  lids  of  the  two  chambers.  The  absence  of  the 
shaft  rotation  ensures  a  small  leakage  of  the  system  which  was 
further  reduced  by  using  an  additional  O-ring  (Seals  DK,  Silicon) 
along  the  shaft  axis.  The  two  chambers  were  filled  with  identical 
aqueous  solutions  by  first  applying  vacuum  (Vacuumbrand  DK  ro¬ 
tary  pump)  in  both  chambers  (valves  V3,  V6  and  V7  close,  all  the 
others  open),  and,  subsequently  by  opening  the  outlet  valves  (VI- 
V5  close,  V6  and  V7  open),  which  were  beforehand  connected  with 
a  solution  reservoir.  After  the  system  reached  the  thermal  stability 
the  low  pressure  chamber  was  connected  to  a  calibrated  capillary 
and  pressure  was  applied  in  the  high  pressure  chamber  by  using  an 
Ar  gas  cylinder  (V3  and  V5  open,  V7  close),  and  was  monitored,  in 
each  cell,  by  the  use  of  pressure  transducers  (Druck  DE,  range  1  — 
11  bar).  The  glass  reactors  were  designed  to  work  above  6  bar  in 
pressure,  and  in  order  to  avoid  any  unexpected  increase  two 
additional  rupture  discs  (Biichi  CFI,  wet  part  in  graphite  and  PFA 
6  bar)  were  connected  to  the  two  cells. 

The  active  area  for  the  membrane  transport  was  3.14  cm2.  To 
avoid  any  deformation  of  the  membrane  due  to  the  applied  pres¬ 
sure,  a  screen  (Sartorius  US,  coarse  grid  Teflon  coated  stainless 
steel)  has  been  used  in  the  low  pressure  compartment.  The 
streaming  potential  was  measured  through  the  use  of  two  Ag/Ag- 
C1  electrodes  (0.25  mm  in  diameter)  placed  close  to  the  membranes 
surfaces  and  mounted  through  the  two  membranes  holders  with 
PEEK  1/16  inch  connector  (Mikrolab  DK)  and  purposely  designed 
Teflon  ferrules.  The  electrodes  were  prepared  from  0.25  mm 
diameter  silver  rod  (Alfa  Aesar  DE,  99.9985%  metal  basis)  following 
the  same  procedure  reported  in  a  previous  publication  [18].  The 
hydraulic  permeability  was  monitored  following  the  movement  of 
the  liquid  meniscus  over  time  in  the  calibrated  capillary  (Teflon, 
inner  diameter  0.75  mm). 

All  the  fitting,  valves  and  tubing  were  in  Teflon  and  bought  from 
Swagelok,  with  the  exception  of  the  GL  connection  of  the  cells  (PFA) 
and  the  vacuum  fittings  and  valves  which  were  purchased  form 
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Fig.  1.  Schematic  of  the  experimental  set-up  used  in  this  work  for  hydraulic  permeability  and  streaming  potential  tests. 


Bola,  DE  and  Edwards,  UK  respectively.  All  the  signals  were  ac¬ 
quired  via  NI-9219  acquisition  card  (National  Instruments)  at  1  Hz 
and  monitored  and  recorded  with  Labview. 

2.3.  Ion  conductivity  through-plane  and  in-plane  apparatuses 

Ion  conductivity  was  measured  in  two  different  set-ups.  In 
particular  through-plane  measurements  were  carried  out  in  the 
same  electrochemical  cell  used  in  a  previous  publication  [18] 
placed  in  an  incubator  (Binder  DE,  maximum  temperature 
100  °C);  the  temperature  of  the  system  was  monitored  with  two 
thermocouples  type  I<  (Omega  Us,  1/32  inch).  Unlike  the  previous 
study,  the  electrodes  were  mounted  in  spring-loaded  holders 
(sketched  schematically  in  Fig.  2)  to  reduce  the  distance  between 
electrodes  and  membrane  surfaces.  In  this  case  to  guarantee  a  more 
planar  surface  the  silver  rod  was  mounted  in  a  tubing  sleeve  (Idex 
US,  0.4  and  0.75  mm  I.D.  and  O.D.  respectively),  glued  and  then 
polished  with  AIO2  polishing  films  of  decreasing  grit  sizes.  After  this 
process  the  exposed  Ag  surface  was  cleaned  with  isopropyl  alcohol 
and  electroplated  following  the  procedure  in  Ref.  [18].  Experiments 
were  performed  with  and  without  screens  (membrane  active  area 
of  71  and  220  mm2  respectively)  resulting  in  a  difference  in  the  ion 
conductivity  of  approximately  6%  (2.53  and  2.68  S  m-1  with  and 
without  screen  respectively  at  25  °C).  The  active  area  for  the  ex¬ 
periments  with  the  screen  was  analyzed  from  microscope  pictures 
with  Matlab.  The  membrane  holders  and  screens  were  3D-printed 
in  PLA  (Makertbot  Replicator  2),  however,  with  this  material  the 
electrodes  holders  could  not  be  used  at  temperature  above  45  °C. 
Alignment  and  distance  between  electrodes  was  checked  by 


measuring  the  resistance  between  the  electrodes  in  two  different 
solutions  (0.03  and  0.12  M  LiCl  solutions)  mounted  in  the  mea¬ 
surement  cell  without  a  membrane.  From  the  measured  resistances 
we  estimated  the  misalignment/distance  between  the  electrodes  to 
be  approximately  2—3  pm. 

The  in-plane  ion  conductivity  measurements  were  carried  out 
by  a  four-point  method  by  placing  the  electrodes  in  the  same  plane 
of  the  membrane  and  clamping  them  between  two  glass  plates.  The 


Fig.  2.  Spring-loaded  holders  for  electrodes  used  in-through-plane  ion  conductivity 
measurements. 
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Table  1 

Nafion  117  properties  measured  in  the  present  work:  membrane  thickness  (5),  hydraulic  permeability  (kh).  streaming  potential  coefficient  (v),  through-plane  (a± )  and  in-plane 
(cr || )  ion  conductivity,  and  the  calculated  figure-of-merit  ((3 Ek).  electrokinetic  efficiency  (?7ei<).  hydraulic  permeability  without  electroviscous  effects  (4)  and  power  density  ( P/A ) 
in  the  temperature  range  25—70  °C  and  in  0.03  M  LiCl  solutions.  (3Ek,  rjEK,  k ^  and  P/A  were  calculated  by  using  the  through-plane  ion  conductivity  values.  Power  densities  have 
been  calculated  considering  a  trans-membrane  pressure  difference  of  10  bar.  Notes:  (a)  interpolated  and  (b)  extrapolated  values  calculated  from  the  activation  energy  of  the 
process. 


T°C 

5  lO"6  m 

kh  1017  m2  s"1  Pa-1 

v  10"9  VPa"1 

<r±  S  m  1 

(T||  S  m  1 

P  EK 

r?Ei<% 

4  1017  m2  s"1  Pa-1 

P/A  W  ITT2 

25 

201 

4.43 

4.1 

2.56 

4.3a 

1.04 

17.68 

8.86 

0.33 

30 

201 

5.28a 

4.4 

3.00a 

4.9a 

1.13 

18.70 

11.25 

0.41 

35 

202 

6.06 

4.7a 

3.48a 

5.4 

1.22 

19.71 

14.18 

0.51 

40 

202 

7.66a 

5.0a 

4.02 

6.2 

1.32 

20.73 

17.77 

0.63 

50 

202 

11.08 

5.9 

5.27b 

7.5 

1.52 

22.75 

27.44 

0.95 

60 

204 

15.87a 

6.2a 

6.81b 

8.2 

1.75 

24.73 

41.48 

1.39 

70 

206 

20.40 

6.7 

8.67b 

9.5 

1.98 

26.68 

61.44 

1.99 

distance  between  the  reference  electrodes  was  12  mm,  and  the 
sample  width  10  mm.  The  sample  holder  was  then  immersed  in 
0.03  M  LiCl  solution  in  a  jacked  borosilicate  reactor  and  heated  at 
the  test  temperature  by  using  a  water  bath  (Eco  Gold  Lauda  DE). 

In  both  cases  the  ion  current  (±1  pA  to  ±2  mA)  was  generated 
with  an  Agilent  U2722A  source  measure  unit  and  the  potential  was 
acquired  with  a  NI-9219  card  and  recorded  with  Labview  Signal- 
Express  at  1.95  Hz. 


Fig.  3.  Upper  panel:  volume  fluxes  as  a  function  of  the  trans-membrane  pressure 
difference  in  the  temperature  range  25-70  °C  in  0.03  M  LiCl  solution.  Lower  panel: 
Arrhenius  plot  of  the  hydraulic  permeability  and  comparison  with  literature  data 
[18,36-39]  for  Nafion  membranes  in  Li+,  Na+  and  H+  form. 


3.  Results 

3.1.  Hydraulic  permeability 

The  upper  panel  of  Fig.  3  shows  the  volume  flux,  jv,  as  a  function 
of  the  applied  pressure  difference  in  the  temperature  range  25- 
70  °C.  Pressure  was  applied  after  the  temperature  of  the  system  was 
stabilized  and  the  temperature  difference  between  the  high  and 
low  pressure  chambers  was  within  0.2  °C.  In  both  half  cells  the 
stirring  rate  was  set  to  the  maximum  value  of  900  rpm  which  at  the 
same  time  gave  better  temperature  stability  as  well  as  diminishing 
the  effect  of  the  concentration  polarization.  As  reported  in  a  pre¬ 
vious  work  [18],  the  presence  of  a  larger  concentration  polarization 
boundary  layer  resulted  in  a  lower  xH,  hence  the  values  reported  in 
this  study  represent  an  upper  limit,  and  should  be  considered  as  a 
good  assessment  of  kh  in  the  real  operating  conditions  in  a  power 
generation  device.  It  should  also  be  noted  that  the  use  of  these 
values  represents  a  conservative  way  to  calculate  the  efficiency  of 
the  system.  The  experiments  were  performed  in  pressure  range 
1  bar  to  approximately  3  bar  and  jv  was  recorded  after  the  steady 
state  flux  was  reached  (approximately  half  hour).  The  /cH  were 
calculated  from  the  slope  of  the  linear  regressions  reported  in  Fig.  3, 
constrained  to  pass  through  the  origin,  and  the  membrane  thick¬ 
ness  in  the  swollen  state.  Membrane  thickness  and  kh  at  several 
temperatures  are  reported  along  with  the  other  measured  prop¬ 
erties  in  Table  1.  To  check  the  consistency  of  the  measured  mem¬ 
brane  thicknesses  an  additional  measurement,  not  reported  in 
Table  1,  was  carried  out  above  the  temperature  range  used  in  the 
present  work  (100  °C)  finding  5  =  211  pm  (swelling  ratio  of  15.2%) 
in  rather  good  agreement  with  the  measurement  of  Sakai  et  al.  [32] 
8  =  220  (swelling  ratio  of  18.9%)  found  at  100  °C  in  Nafion  in  H+ 
form.  Nafion  in  H+  and  Li+  forms  have  similar  water  uptake  [38] 
and  consequently  should  have  similar  dimensional  swelling.  The 
value  for  kh  found  at  25  °C  of  4.41  10  17  m2  s-1  Pa-1  is  somewhat 
larger  with  respect  to  the  one  found  at  22  °C  in  the  previous  work 
3.1 10-17  m2  s_1  Pa-1  [18].  However,  the  slightly  higher  temperature 
in  the  present  study  can  explain  about  half  the  difference.  We 
believe  that  the  rest,  which  is  of  the  same  order  as  the  resolution  of 
the  experiment,  could  be  explained  by  a  better  stirring  in  the 
present  experimental  setup  and,  as  a  consequence,  a  thinner 
boundary  layer. 

In  the  lower  panel  of  Fig.  3,  kh  is  reported  as  a  function  of 
temperature  in  an  Arrhenius  plot  and  compared  with  literature 
data  for  Nafion  (Nafion  117  and  Nafion  115)  in  Li+  [18],  Na+  [18]  and 
H+  [36-39]  form.  Data  at  25  °C  collected  by  Evans  et  al.  [38]  in  H+ 
form  are  reported  to  highlight  the  difference,  almost  four-fold, 
between  xH  for  the  treated  and  untreated  Nafion  samples.  Hence 
the  difference  in  the  hydraulic  permeability  found  between  our 
previous  and  current  study  can  be  also  partially  attributed  to  the 
slight  difference  in  the  pretreatment  protocol  followed. 
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The  apparent  activation  energy  for  the  hydraulic  permeability  of 
Nation  in  0.03  M  LiCl  solutions,  EKh,  was  29.4  kj  mol-1,  and  to  the 
best  of  our  knowledge  it  has  not  been  reported  by  any  other  groups 
in  the  literature.  However,  the  value  found  in  the  present  study  is 
within  the  one  calculated  by  Duan  et  al.  [37],  24.9  kj  mol-1  (T  =  23- 
80  °C),  and  the  one  of  Adachi  et  al.  [36]  of  37.8  kj  mol-1  (T  =  21- 
70  °C)  for  Nation  in  H+  using  pure  water. 


3.2.  Streaming  potential 

As  reported  in  one  of  our  previous  studies  [18]  measurements  of 
the  streaming  potential  coefficient  with  and  without  stirring  gave 
the  same  values  within  5%.  In  the  present  study  we  use  the  method 
without  stirring  only:  the  advantage  of  using  this  method  is  a  more 
stable  signal  of  the  electrodes,  and,  since  the  experiments  were 
conducted  for  a  maximum  of  200  s,  no  significant  change  in  the 
temperature  of  the  system  was  observed.  With  this  method  a 
concentration  polarization  layer  is  built  up  and  contributes  to  the 
streaming  potential.  However,  the  intrinsic  membrane  contribution 
is  time  independent,  whereas  the  contribution  from  concentration 
polarization  has  yjt  behavior  being  the  diffusion  of  ions  in  the  so¬ 
lution  a  diffusion-limited  process.  Hence,  the  streaming  potential 
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Fig.  4.  Upper  panel:  potential  difference  as  a  function  of  the  trans-membrane  pressure 
difference  in  the  temperature  range  25-70  °C  in  0.03  M  LiCl  solution.  Lower  panel: 
Arrhenius  plot  of  the  streaming  potential  coefficient  and  comparison  with  literature 
data  for  Nafion  in  0.03  M  LiCl  solution  at  ambient  temperature  [18,  41]. 


coefficient  can  be  found  from  fitting  A </>(t)  =  vAp  +  AJt  to  the  data 

[40]. 

The  upper  panel  of  Fig.  4  shows  the  measured  streaming  po¬ 
tential  in  the  temperature  range  25-70  °C  as  a  function  of  the 
applied  trans-membrane  pressure  difference.  Stirring  is  applied 
until  the  measurement  temperature  is  reached,  but  ceased  before 
the  measurement  is  started.  From  the  linear  regression,  constrained 
to  pass  through  the  origin,  of  the  streaming  potential  coefficients,  v 
are  found  at  four  different  temperatures  and  they  are  shown  in  an 
Arrhenius  plot  in  the  lower  panel  of  Fig.  4.  v  varied  from 
4.12  1(T9  V  Pa”1  to  6.7  1(T9  V  Pa-1  at  25  °C  and  70  °C.  respectively. 
The  apparent  activation  energy  for  v  (Ev)  was  approximately 

9.3  kj  mol-1.  By  considering  this  value  for  the  activation  energy  it 
was  possible  to  extrapolate  v  to  lower  temperatures.  Specifically  at 
22  °C  the  extrapolated  value  of  4.04 10-9  V  Pa-1  was  fully  consistent 
with  our  previous  work  [18]  (6%  higher).  Compared  to  the  value  of 
2.95  10-9  V  Pa-1  found  by  Xie  and  Okada  [41]  the  difference  is 
larger  and  might  be  due  to  a  membrane-batch  dependence  and/or 
differences  in  the  pretreatment  of  the  membrane  as  already  dis¬ 
cussed  in  Ref.  [18]. 

A  direct  comparison  of  Ev  with  literature  data  is  not  possible 
because  of  lack  of  experimental  data.  Nonetheless,  streaming  po¬ 
tential  experiments  were  performed  by  Okada  et  al.  [42]  in  Nafion 
117  by  using  0.03  M  NaCl  and  0.03  M  HC1  in  the  temperature  range 
25-60  °C  and  25-80  °C  for  the  two  solutions,  respectively.  The 
calculated  Ev  for  Nafion  in  Na+  and  H+  form  were  3.2  and 

7.3  kj  mol-1  respectively  in  the  temperature  range  25-60  °C.  As 
could  have  been  anticipated  the  latter  value  is  comparable  with  the 
one  retrieved  in  the  present  study,  since  Nafion  in  H+  and  Li+  forms 
have  similar  water  uptake.  However,  the  streaming  potential  co¬ 
efficient  measured  by  Okada  et  al.  42]  in  Nafion  in  H+  form  at  80  °C 
was  lower  than  the  one  at  60  °C.  If  also  this  data  point  is  included 
in  the  regression  for  the  activation  energy  calculation 
Ev  ~  3.6  kj  mol-1,  a  value  2  Vi  times  lower  than  the  one  found  in  the 
present  work.  No  explanation  for  the  anomalous  behavior  of  v  at 
80  °C  was  reported  by  Okada  et  al.  [42].  As  mentioned  previously 
the  difference  with  respect  to  the  data  in  Ref.  [42]  might  be 
attributed  to  a  membrane-batch  dependence  or  a  different  mem¬ 
brane  pretreatment.  This  interpretation  is  supported  by  the  low 
value  of  the  activation  energy  for  the  water  permeability, 
11.6  kj  mol-1,  in  Nafion  H+  found  by  Okada  et  al.  [42]. 

3.3.  Ion  conductivity 

Results  for  ion  conductivity  measurements  are  reported  in  Fig.  5. 
The  upper  panel  shows  the  ion  conductivity  (a±)  measured 
through-plane  with  a  four-electrode  experimental  setup  in  the 
temperature  range  22-40  °C.  The  applied  current  was 

between  -1  mA  and  1  mA.  Measurements  at  22  °C  were  carried  out 
to  have  a  comparison  with  data  reported  by  Kilsgaard  et  al.  [18]  in 
the  same  membrane  and  solution.  In  this  case  the  agreement  is 
excellent  (2.32  S  m-1  in  the  present  study  and  2.36  S  m-1  in 
Ref.  [18]).  As  reported  in  a  previous  study  [18]  literature  values  of  a 
measured  in  LiCl  solutions  ranged  between  0.6  S  m-1  and 
approximately  3  S  m-1  [43-47]  and  the  rather  marked  differences 
can  be  explained  by  the  method  used,  the  different  water  uptakes 
of  the  membrane  and/or  membrane  pretreatments. 

As  stated  earlier,  measurements  of  a±  are  limited  to  tempera¬ 
tures  below  45  °C  in  which  we  calculate  an  activation  energy  (Ea± ) 
of  22.9  kj  mol-1.  This  value  is  close  to  23.5  kj  mol-1  measured  by 
Okada  et  al.  [42]  in  Nafion  Na+  form  with  AC  impedance  method. 

Measurements  of  the  in-plane  conductivity  (cr|| )  were  done  in 
order  to  reach  higher  measurement  temperatures.  Data  are  shown 
in  the  lower  panel  of  Fig.  5.  Due  to  the  different  membrane  sample 
geometry  used  in  the  in-plane  experiments,  the  applied  current 


J.  Catalano,  A.  Bentien  /  Journal  of  Power  Sources  262  (2014)  192-200 


197 


Fig.  5.  Ion  conductivity  measurements  in  0.03  M  LiCl  solutions.  Upper  panel:  4- 
electrode  through-plane  tests  (22-40  °C).  Lower  panel:  4-electrode  in-plane  tests 
(22-70  °C). 


an  index  of  the  degree  of  anisotropy  of  the  membrane:  in  the 
present  case  (7||/t7±  1 22oC  =  1.7  while  values  of  3.6,  2.5  and  1.2  were 
found  by  Gardner  and  Anantaraman  [49,50],  Eldab  et  al.  [48]  and 
Soboleva  et  al.  [51],  respectively.  As  seen  in  Fig.  6  the  degree  of 
anisotropy  decreases  as  the  temperature  of  the  system  is  increased. 
This  effect  can  be  explained  because  by  increasing  the  temperature 
the  water  fraction  inside  the  polymer  increases  reducing  the  in¬ 
fluences  of  a  preferential  orientation.  It  should  also  be  noted  that 
isotropic  conductivity  (e.g.  Ref.  [53  )  and  reverse  anisotropy  [54]  (in 
Nation  112)  were  also  found  in  extruded  Nation  membranes. 

The  difference  in  Effjj  andF^  might  be  a  consequence  of  the 
anisotropy  of  the  membrane  considered,  hence  a±  (as  well  as  its 
activation  energy)  should  be  used  in  the  calculation  of  the  effi¬ 
ciency  for  the  electrokinetic  power  generation  since  it  was 
measured  in  the  same  direction  of  the  ion-transport.  Therefore  the 
extrapolated  values  of  a±  in  the  temperature  range  50-70  °C  will 
be  used  for  the  figure-of-merit  calculations.  It  should  also  be  noted 
that  the  use  of  a  ±  in  the  calculation  of  the  figure-of-merit  repre¬ 
sents  a  conservative  way  for  evaluating  the  efficiency  of  the  system. 

4.  Discussion 

In  Table  1  all  the  properties  measured  in  the  present  work  are 
reported,  namely  5 ,  /cH,  v,  a ±  and  o\\  and  the  calculated  figure-of- 
merit ,  electrokinetic  energy  conversion  efficiency,  non- 
electroviscous  hydraulic  permeability  (k*h)  [16]  and  power  density 
(P/A)  in  the  temperature  range  25—70  °C.  When  the  property  was 
not  measured  at  the  specific  temperature,  the  property  has  been 
calculated  by  interpolation  of  measured  data.  As  reported  previ¬ 
ously,  the  only  exception  is  the  through-plane  ion  conductivity,  for 
which  at  temperature  equal  or  above  50  °C  the  value  has  been 
calculated  by  extrapolation.  The  calculated  efficiency  and  power 
density  are  also  reported  as  a  function  of  the  temperature  in  Fig.  7.  It 
is  worth  to  note,  that  the  efficiency  at  room  temperature,  17.7%,  is 
close  to  the  one  found  in  one  of  our  previous  work  (20  ±  2%)  [18]  at 
22  °C.  At  70  °C,  the  calculated  efficiency  reaches  the  rather  high 
value  of  26.7%,  approximately  50%  higher  than  that  at  25  °C  with  a 
power  density  of  ca.  2  W  m-2.  By  using  the  calculated  activation 


range  was  greatly  reduced  and  it  was  in  the  range  -20  pA— 20  pA. 
For  cr ||  we  calculate  an  activation  energy  £ff||  =  15.4  kj  mol-1  as  seen 
in  Fig.  6  that  is  somewhat  lower  than  Ea± .  In  addition  it  is  seen  that 
c ||  is  systematically  higher  with  respect  to  (e.g.  at  22  °C 
o ||  =  3.8  S  m-1  while  a±  =  2.32  S  m-1)  hence  o\\  represents  an 
upper  limit  of  the  cation  conductivity  of  the  membrane.  Proton 
conductivity  experiments  performed  along  and  perpendicularly  to 
the  membrane  extrusion  axis  reported  by  several  authors  [48-51  ] 
showed  an  anisotropic  proton  conductivity  of  Nation  which  has 
been  ascribed  either  to  the  preferential  orientation  at  the  mem¬ 
brane  interface  or  to  the  orientation  of  Nation  chains  along  the 
extrusion  axis.  It  could  be  suspected  that  the  lower  value  of  a ± 
could  be  a  systematic  error  from  the  experimental  setup  where 
resistive  contributions  from  the  electrolyte  solution  would 
decrease  the  apparentcr±.  Nonetheless,  we  have  confidence  that 
the  difference  between  a y  and  a±  is  intrinsic.  Measurements  of  a± 
are  reproducible  (not  shown)  and  as  mentioned  earlier,  measure¬ 
ments  in  the  cell  without  membrane  show  that  misalignment  of 
the  electrodes  is  approximately  2-3  pm,  whereby  electrolyte 
contribution  to  the  total  cell  resistance  becomes  negligible. 

Experiments  made  in  Nation  membrane  without  any  preferen¬ 
tial  orientation  [52]  (casted  Nation  NR-212)  reporting  isotropic 
conductivity  substantiate  the  idea  of  the  anisotropy  due  to  the 
extrusion  process.  Also  the  ratio  between  g\\/gj_  was  proposed  as 
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Fig.  6.  Arrhenius  plot  of  the  ion  conductivity  measurements  performed  in  both  con¬ 
figurations  (through-plane  and  in-plane)  in  0.03  M  LiCl  solutions.  For  comparison  data 
from  Kilsgaard  et  al.  [18]  retrieved  in  Nation  117  and  same  solutions  have  been  re¬ 
ported  (through-plane).  Data  from  Elabd  et  al.  [48]  for  proton  conductivity  in  Nation 
117  measured  with  through-plane  and  in-plane  configurations  are  also  shown. 
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Fig.  7.  Calculated  electrokinetic  efficiency  and  power  density  of  Nation  117  membrane 
in  0.03  M  LiCl  solution  as  a  function  of  the  temperature.  Power  densities  have  been 
calculated  considering  a  trans-membrane  pressure  difference  of  10  bar.  Lines  are  a 
guide  for  the  eyes. 


energies,  reported  in  Table  2  along  with  the  pre-exponential  fac¬ 
tors,  for  the  processes  involved  one  can  write: 


ftx(r) 


Vq-OQ 

kh,o 

E'o'gQ 

kH,0 


exp 

exp 


2  Ev  +  Ea±  —  EKh\ 

RT  J 

RT  J 


(1) 


In  which  represents  the  apparent  activation  energy  of  the 
electrokinetic  process  and  in  the  case  of  Nation  is  approximately 
12.2  kj  mol"1. 

The  power  density  PI  A  was  calculated  as  the  work  produced  per 
unit  time  and  area  by  the  liquid  in  the  membrane  assuming  a 
transmembrane  pressure  difference  of  10  bar  using  the  formula 
P/A  =  KH-Ap2/<5( l+2jtfEI<);  the  latter  expression  is  straightfor¬ 
ward  derived  from  the  phenomenological  equations  and  does  also 
include  electro-viscous  effects  (V^ek)  which  increase  the  apparent 
hydraulic  permeability  and  thereby  power  density,  during  opera¬ 
tion.  Still  the  calculated  power  density  of  electrokinetic  power 
generation  is  low  compared  to  fuel  cells  [55].  However,  it  can  most 
likely  be  increased  up  to  more  than  100  W  m-2  by  using  higher 
pressure  difference,  thinner  or  supported  composite  membranes 
with  a  very  thin  selective  ion  conductive  layer.  Finally  Table  1  also 
report  the  non-electroviscous  hydraulic  permeability  [16], 
%  =  kh(1  +fe<)»  which  physically  represents  the  hydraulic 
permeability  when  a  current  density  is  allowed  to  circulate  through 
the  membrane.  The  electroviscous  effect  arises  because  of  the 
presence  of  charges  imbedded  in  the  membrane  which  influence 
the  permeation  of  the  liquid  in  the  water  channels;  if  the  charges 
are  free  to  move  the  electroviscous  effect  would  not  be  present  and 
consequently  the  membrane  permeability  would  increase. 


Table  2 

Pre-exponential  factors  and  activation  energies  for  the  transport  properties 
measured  in  this  work  between  ambient  temperature  and  70  °C.  Notes:  (a) 
measured  between  22  and  40  °C. 


Pre-exponent 

Activation  energy  kj  mol  1 

kh 

6.22  10-12 

m2  s"1  Pa"1 

29.4 

V 

1.80  10"7 

V  Pa-1 

9.3 

2.70  104 

S  ITT1 

22.9 

ff\\ 

2.10  103 

S  m  1 

15.4 

0EK 

1.40  102 

12.2 

The  electrokinetic  transport  properties  can  be  described,  in  a 
simplified  way,  in  terms  of  a  parallel  array  of  straight  channels,  see 
e.g.  Refs.  [28,29]  for  a  detailed  derivation  of  the  equations: 


<?wr2# 

8  V 


(2) 


1  ftf 

F8t] 


(3) 


f2xo 


(4) 


where  <pw,r,  $,  ??,/,  Fand  X  are  the  volumetric  fraction  of  water  in  the 
polymer,  average  radius  of  the  pore,  the  tortuosity  of  the  matrix,  the 
viscosity  of  the  aqueous  solution,  the  frictional  coefficient,  the 
Faraday  constant  and  the  charge  density  per  unit  volume  of  the 
membrane,  respectively. 

The  /ch  is  estimated  by  the  average  radius  dimension  in  which 
the  volume  flux  takes  place,  not  considering  the  friction  between 
the  water  passing  through  the  membrane  and  the  counter-ion  that 
are  considered  fixed.  It  also  worth  to  note  that  in  general  the  ion 
conductivity  can  be  expressed  as  the  sum  of  the  terms  related  with 
the  bulk  concentration  of  ions  in  the  solution  and  the  counterions 
shielding  the  surface  charge.  However  in  the  presence  of  very  dilute 
solutions  the  contribution  of  the  former  term  to  the  total  ion  con¬ 
ductivity  can  be  considered  negligible  [56]. 

Eqs.  (2)  and  (3)  can  be  separately  used  to  have  two  independent 
estimations  of  the  pore  size  of  a  charged  membrane  even  if  some 
additional  information  on  the  membrane  characteristics  are 
needed:  i.e.  volumetric  fraction  of  water  and  the  tortuosity  of  the 
water  channels  in  the  polymeric  matrix.  However,  in  first  approx¬ 
imation,  $  can  be  set  equal  to  one  as  in  Ref.  [57]  and  the  volumetric 
fraction  of  water  can  be  derived  by  the  knowledge  of  the  linear 
swelling  under  the  hypothesis  of  isotropic  volume  dilation.  By 
considering  the  equation  for  the  frictional  coefficient,  /  =  RT/Dy  in 
which  Dy  =  Du,oexp(-Fdiff  ulRT)  represents  the  diffusion  coefficient 
of  Li+,  it  is  thus  possible  to  estimate  the  average  pore  size  of  the 
parallel  array  using,  as  further  simplification,  for  r/  the  value  in  pure 
water  and  for  Dy.o  and  Edifr  u  the  one  in  dilute  aqueous  solutions 
2.76  10-6  m2  s-1  and  19.5  kj  mol-1,  respectively  [58-60]. 

These  calculations  have  been  reported  in  terms  of  average  pore 
size  of  the  water  channel  in  Fig.  8  as  a  function  of  the  temperature. 
The  pore  size  diameter  increases  from  roughly  2.0  to  2.8  nm  in  the 
temperature  range  25-70  °C  assuming  similar  values  for  the  two 
independent  measurements.  This  latter  evidence  substantiates  the 
internal  consistency  of  the  hydraulic  permeability  and  streaming 
potential  coefficient  values  found  in  the  present  work.  Indeed  by 
combining  Eqs.  (2)  and  (3): 


kr  kh,o 
— = —  exp 
v  v0 


RT 


Dupexp 


RT  J 


(5) 


it  is  apparent  that  the  temperature  dependence  of  the  ratio  be¬ 
tween  the  hydraulic  permeability  and  streaming  potential  coeffi¬ 
cient  should  be  similar  to  the  activation  energy  of  the  diffusion 
coefficient  of  lithium  ion  (neglecting  the  slight  temperature 
dependence  of  <pw).  From  the  experiments  reported  in  the  present 
work  EKh  —  Ev  =  20.1  kj  mol  1  and  resulted  similar  to  Edifr 
u  =  19.5  kj  mol-1. 

It  should  be  pointed  out  that  the  calculated  pore  dimensions  are 
not  representative  of  the  real  structure  of  Nafion  membranes  but 
can  be  used  to  compare  the  performances  of  different  charged 
membranes.  The  same  figure  also  reports  results  from  Schmidt- 
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Fig.  8.  Esteemed  average  pore  diameters  from  the  analysis  of  the  streaming  potential 
and  hydraulic  permeability  measurements  assuming  a  parallel  array  of  straight 
channels.  For  comparison  sake  are  also  reported  the  calculated  size  from  analysis  of 
SAXS  experiments  in  Nation  H+  form  at  20%  vol  of  water  by  Schmidt-Rohr  and  Chen 
[27]  and  from  hydraulic  permeability  in  pure  water  by  Duan  et  al.  [37]  in  Nation  115. 


Rohr  and  Chen  [27]  retrieved  from  SAXS  data  on  Nation  membrane 
by  using  parallel  water-channel  model  simulations.  In  Ref.  [27]  the 
authors  used  long  straight  parallel  channels  in  cylindrical  inverted 
micelles  model,  able  to  reproduce  the  characteristic  SAXS  features 
at  different  hydration  levels.  In  that  study  channels  diameters  from 
1.8  to  3.5  nm,  with  an  average  of  2.4  nm,  were  found  at  20%  vol. 
(A  =  6.7molH2o /molsosH  )  so  at  a  hydration  level  well  below  the  one 
used  in  this  work.  However  the  pore  diameters  calculated  in  the 
present  work  are  of  the  same  order  of  magnitude,  compatible  with 
the  believed  structure  of  Nation  [22-24]  and  with  the  calculated 
average  dimensions  from  hydraulic  permeability  measurements  by 
Duan  et  al.  [37]. 

In  addition  Eq.  (4)  can  be  used  to  critically  discuss  the  ion 
conductivity  values  used  in  the  electrokinetic  efficiency  calcula¬ 
tions  since  it  only  requires  the  knowledge  of  the  diffusion  coeffi¬ 
cient  of  Li+  and  the  charge  density  per  unit  volume  of  the  swollen 
membrane.  This  latter  quantity  can  be  expressed  as  a  function  of 
the  volume  dilation,  the  ion  exchange  capacity,  /EC,  and  the  dry 
polymer  density  ppoi:  X  =  IEC  x  ppoi/(l  +  AV/V).  The  calculated 
values  of  a  are  higher  than  the  one  extrapolated  from  the  experi¬ 
mental  results  (48.3%  and  17.3%  higher  at  50  and  70  °C  respectively) 
hence  the  latter  represent  a  conservative  way  for  evaluating  the 
efficiency  of  the  system. 

Finally  the  set  of  equations  can  be  substituted  in  Eq.  (1) 
obtaining: 


fr 2  X 


(6) 


This  equation  can  be  used  as  a  first  approximation  of  the  elec¬ 
trokinetic  efficiency  of  an  ionomeric  membrane  by  the  knowledge 
of  its  pore  size,  dry  density,  volume  dilation  and  ion  exchange  ca¬ 
pacity.  Also  it  substantiates  the  idea  that  charged  membranes  with 
larger  pore  size,  such  as  collodion  cation  exchange  membranes  or 
modified  Nafion,  would  have  higher  efficiency  as  well  as  power 
density. 


5.  Conclusions 


The  main  goal  of  the  present  study  was  to  measure  the  elec¬ 
trokinetic  power  generation  efficiency  of  Nafion  membranes  at 
different  temperatures.  In  the  temperature  range  25  °C-70  °C,  the 


activation  energies,  of  three  electrokinetic  transport  properties, 
namely  hydraulic  permeability,  streaming  potential  coefficient  and 
ion  conductivity  in  0.03  M  LiCl  solutions  were  measured 
(29.4  kj  mol-1,  9.3  kj  mol-1  and  22.9  kj  mol-1  respectively).  The 
efficiency  was  evaluated  from  the  electrokinetic  figure-of-merit  and 
showed  a  monotonous  increase  between  about  17.7  and  26.6%  by 
increasing  the  temperature  of  the  system  from  25°Cto70°C.  From 
a  simple  model,  that  considers  the  Nafion  membrane  as  an  array  of 
cylindrical  parallel  channels  in  which  the  transport  takes  place,  it 
was  possible  to  calculate  the  average  channel  diameters.  This  latter 
property  was  shown  to  vary  between  approximately  2.0-2.8  nm  (at 
25  °C  and  70  °C  respectively)  in  rather  good  agreement  with  the 
believed  microstructure  of  Nafion  membrane  in  swollen  state. 
These  calculations  additionally  showed  that  using  membranes  with 
larger  pore  size  could  eventually  result  in  higher  power  generation 
efficiency. 
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Nomenclature 

Du  diffusion  coefficient  of  Li+,  m2  s-1 
Du,0  pre-exponential  factor  for  the  diffusion  coefficient  of  Li+, 
m2  s-1 

E  activation  energy,  J  mol-1 

F  Faraday  constant,  96,485  C  mol-1 

/  frictional  coefficient,  J  s  mol-1  m-2 

I  ion  current,  A 

IEC  ion  exchange  capacity,  mol  kg-1 

jv  volume  flux,  m  s-1 

A p  trans-membrane  differential  pressure,  Pa 

P/A  power  density,  W  m-2 

R  ideal  gas  constant,  J  mol-1  K-1 

r  pore  radius,  m 

T  temperature,  K 

X  charge  density  per  unit  volume,  mol  m-3 

AV/V  volume  dilation 

Greek  symbol 

0 ek  electrokinetic  figure-of-merit 

d  membrane  thickness,  m 

7]  viscosity,  Pa  s 

tjek  electrokinetic  efficiency 

&  tortuosity  of  the  matrix 

K\\  hydraulic  permeability,  m2  s-1  Pa-1 
kr,o  pre-exponential  factor  hydraulic  permeability, 

m2  s-1  Pa-1 

A  solution  uptake,  molH2o(molso3H-)-1 

v  streaming  potential  coefficient,  V  Pa-1 

vo  pre-exponential  factor  streaming  potential  coefficient, 

V  Pa-1 

p  density,  kg  m-3 

a  ion  conductivity,  S  m-1 

Go  pre-exponential  factor  ion  conductivity,  S  m-1 

Gq  charge  density,  C  m-2 

cp  volumetric  fraction 

</>  streaming  potential,  V 
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Superscripts 

*  referred  to  intrinsic 

Subscripts 

diff  Li  referred  to  diffusion  coefficient  of  Li+ 

pol  referred  to  polymer  in  dry  condition 

w  referred  to  water 

/?EK  referred  to  electrokinetic  figure -of-merit 

kh  referred  to  hydraulic  permeability 

v  referred  to  streaming  potential  coefficient 

a  referred  to  ion  conductivity 

J_  referred  to  through-plane  experiments 

||  referred  to  in-plane  experiments 
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